Scattering laser light on cold atoms: multiple scattering signals from single-atom 

responses 



Thomas Wellens,"'^ Tobias Geiger/ Vyacheslav Shatokhin,'^' ^ and Andreas Buchleitner-"^ 

^Institut fur Physik, Albert-Ludwigs-Universitdt Freiburg, 
Hermann- Her der-Str. 3, D-79104 Freiburg, Germany 
^B. I. Stepanov Institute of Physics NASB, 220072 Minsk, Belarus 
(Dated: September 30, 2009) 

We deduce the coherent backscattering signal from two distant laser-driven atoms using single- 
atom equations. In contrast to the standard master equation treatment, this new approach is 
suitable for the generalization to a large number of atomic scatterers. 
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Cooled atomic gases have been successfully used in 
experimental studies of coherent backscattering (CBS) 
of laser light [1 , and are considered to be a promis- 
ing medium for the observation of Anderson localization 
[2] with photons. Apart from their fundamental inter- 
est, these interference phenomena in disordered systems 
are important for practical applications such as random 
lasers 0] and quantum memory for light [4J. Inasmuch 
as multiple scattering of photons is crucial for the emer- 
gence of weak and strong localization [S] in disordered 
atomic clouds, while it can play a detrimental role for 
the storage of quantum information [6j, an appropriate 
theory taking into account all relevant effects of atom- 
light interactions is in order. 

The theory of multiple scattering in dilute media that 
consist of a disordered collection of discrete scatterers 
relies on the division of the total scattering process into 
single scattering events. These are described by a scat- 
tering operator which is usually assumed to be identi- 
cal for all the individual scatterers, and thereby serves 
as the fundamental building block for the total multiple 
scattering process [7]- This approach, however, fails in 
the case of atomic scatterers, for the following two rea- 
sons: first, due to the saturation of the atomic transition 
with increasing intensity of the incident light, the atoms 
scatter Hght nonUnearly. The outgoing field (i?out) is 
not proportional to the incoming one (ii^in), and there- 
fore cannot be described by a linear scattering operator 
of the form Eont — TEi^. Second, the hght scattered 
by near-resonant atoms exhibits fluctuations due to the 
quantum mechanical coupling of the atoms to the electro- 
magnetic vacuum. This means that even if the incident 
laser light i?i„ is perfectly coherent, this is not the case 
for the scattered field i?out- In particular, its average in- 
tensity differs from the square of the average field, i.e., 
(l^^outP) > K-Eout)!^- The difference defines the incoher- 
ent (or inelastic) component of the resonance fiuorescence 
intensity. 

So far, no satisfying theory exists for incorporating 
both these effects into a multiple scattering approach. 
Recently [8, 9J, a theory for coherent backscattering by 
nonlinear, classical scatterers was presented, but does not 
take into account any quantum fluctuations due to in- 



elastic scattering. A perturbative approach based on the 
scattering matrix of two photons was proposed in |10j . 
but is only valid if incident light intensity and optical 
thickness of the atomic medium are small. On the other 
hand, standard tools of quantum optics (master equa- 
tions, optical Bloch equations . . . ) are well adapted to 
describe the atom-field interaction for arbitrary intensi- 
ties of the incident field, but are restricted to a small 
number of atoms coupled to each other by photon ex- 
change [TTJ [T^ . This is due to the fact that the dimen- 
sion of the atomic Hilbert space grows exponentially with 
the number of atoms. 

In the present Letter, we show that - for the case of 
two randomly placed atoms with mean distance much 
larger than the laser wavelength - the two perspectives 
of nonlinear multiple scattering and open system dynam- 
ics can be unified within a single approach. Our results 
agree perfectly well with the master equation and quan- 
tum Langevin treatments presented in [T31 Ej . In con- 
trast to [Ol [H], however, the present approach can be 
generalized for a large number of atoms. 

Imagine a cloud of atoms in free space which is ex- 
cited by a coherent laser beam. We are interested in 
the average intensity radiated into different directions 9, 
where the average is taken over the random positions of 
the atoms. Around the backscattering direction 9 — 0, a, 
coherent backscattering (CBS) peak is observed [T], origi- 
nating from the constructive interference of counterprop- 
agating, multiply scattered waves. In the case of a dilute 
and optically thin cloud, the leading contribution to the 
interference signal comes from double scattering events, 
where the laser interacts only with pairs of atoms. (Sin- 
gle scattering does not contribute to CBS, since direct 
and reversed paths are identical in this case.) Hereafter, 
we therefore derive the double scattering CBS contribu- 
tion from a simplified model composed of only two atoms, 
considered stationary for the time for one scattering cy- 
cle to take place, and separated by a distance which is 
large compared to the laser wavelength [TTJ [T^] . For sim- 
plicity, we further consider two-level atoms coupled to 
a scalar photon field, such as to avoid the additional 
(though straightforward) technical overhead needed to 
account for photon polarization and degenerate atomic 
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Figure 1: Double scattering contributions to coherent 
backscattering. The thick arrows represent the intense laser 
beam, and the thin solid and dashed arrows the Fourier com- 
ponents of the scattered fields, E^'^\t) and E^~\t), respec- 
tively, (a) One of the atoms emits a photon with frequency u 
which is scattered by the other atom before it reaches the de- 
tector. This process yields the background intensity Il(i-od) 
of the photodetection signal, (b) One of the atoms emits a 
photon with frequency to, the other a photon with frequency 
uj' = lol -\- uiD — OJ (see text). This process gives rise to the 
coherent backscattering interference peak Icij^o)- 



transitions. 

The situation we consider is summarized in Fig.[l] Two 
atoms scatter photons injected through an intense laser 
beam at frequency lj^ into distinct modes with frequen- 
cies Lo, uj' , and ojd- While cjd is the frequency of the 
detected photons, lo and oj' are the frequencies of the 
(unobserved) photons exchanged between the atoms ~ 
thus mediating the double scattering process. Fig. [TJi de- 
fines the weakly angle-dependent background of the pho- 
todetection intensity, also called "ladder intensity" II in 
the following, whereas Fig. [Tja determines the height Ic 
( "crossed intensity" ) of the coherent backscattering peak 
around = 0, resulting from interference between pho- 
tons travelling along reversed scattering paths. Both II 
and Ic are independent of the positions of the two atoms, 
and thus define, for the case of large mean distance be- 
tween the atoms, the configuration-averaged photodetec- 
tion intensity [iTl [12] . The assumption of large distance 
is required in order to neglect processes where more than 
one photon is exchanged between the atoms. 

The key ingredient for our subsequent treatment is the 
observation that the detected double scattering signal 
can be deduced from the superposition of the fields scat- 
tered by both individual atoms into the detector, each 
atom being subject to a bichromatic driving field £{t) 
composed of the injected laser and the photons scattered 
by the other atom. 



Sit) = Slc 



-iujL i 



+ £le'^^' + v^-^e'^\ (1) 



The classical field £{t) is decomposed into its posi- 
tive and negative frequency components with complex 
amplitudes £l and v^'^\— w^~^*), since the final fig- 
ure of merit, the average stationary intensity (/) — 
limt^ao {E^~Ht)E^~^\t)) measured by the detector, is 
conveniently expressed in terms of positive and negative 
frequency components of the signal radiated by the two 
atoms. 



The single atom response to the bichromatic drive is 
characterized by the associated cross correlation function 
in frequency space, C{uji,lu2) = {E'^^^uji)E'^^^llJ2)) , 
which correlates the, in general, different frequencies of 
the negative (thin dashed lines in Fig. [T|) and positive 
(thin continuous lines in Fig. [TJ frequency components 
^^^LUi) and i?(+)(ti>2) of the emitted field (left atom: 
LUi — and uj2 — uj, right atom: uji = u)' and ^2 — ujd, 
in Fig. hp). Since the emitted fields iJ^^-' and £''^+\ in 
turn, are proportional to the atomic dipole raising and 
lowering operators cr_|_ and (t_ , the frequency correlation 
function can be written as: 



C(wi,W2) 



27r 



(2) 

and is thus given by the solutions cr+ {t) and cr_ {t) of the 
single atom optical Bloch equations under the bichro- 
matic drive, Eq. ([T]). 

The total detected double scattering intensity at fre- 
quency lod is finally given by an integral over the product 
of the cross correlation functions of both atoms, where 
the integral runs over the frequencies of the exchanged 
photons. For an explicit evaluation of these integrals, we 
expand C{uji,lj2) to lowest order in v (+' and This 
perturbative treatment is actually already implied by our 
ansatz of classical scattered fields in Eq. ([T|), which ne- 
glects their quantum statistical properties - encoded in 
the second order correlation functions, i.e. in higher or- 
ders of w^*^ - and is certainly justified in the presently 
considered limit of a dilute gas with large interatomic dis- 
tances and, thus, <C \£l\- Solving Eq. ^ for the 
case of bichromatic driving, Eq. ([T]), and keeping only 
terms up to first order in w*^"*"^ and/or we obtain the 
following four terms: 



dC 



dC 



(5(u;i - W2)fo(wi), (3a) 

5{UJ2 ~ UJi- UJ + LLlL)Pi{uj;UJ2),{'i^) 

5{uii - uo2~ UJ + wl)Pi (w; wi),(3c) 



5{lOi - W2)P2(w; LOi). 



(3d) 



The (5- functions in Eqs. ( 3apd I originate from integrating 
over i+ :— ti+t2 in Eq. (2^ and are thus a consequence of 
time translation invariance or energy conservation. They 
ensure that the negative and positive frequency compo- 
nents of the field are radiated at the same frequency uj 
in Fig. [l^, whereas uj' = uj^ + ujd — uj in Fig. [T]d. 

Po{uj) in Eq. (3a) denotes the usual resonance fiuo- 



rescence spectrum for monochromatic driving including 
both the elastic and inelastic components [16 , and is as- 
sociated with the emission of a fluorescence photon of fre- 
quency UJ by the left atom in Fig.jlji. P2{uj; ujd), defined 
by Eq. ( 3d I , describes the subsequent scattering cross sec- 



tion of this fiuorescence photon by the right atom, into 
the finally detected frequency mode ujo- Consequently, 
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the final expression for the ladder contribution (Fig. [TJi) 
to the total scattered intensity at frequency luo reads 



Il{^d) = \9? I dujP^{uj)P2{uj;uJD) , 



(4) 



where g — exp{iuj]^r) /{ki^r), with kj^ = ujl/c, describes 
the effective coupling strength between the two atoms. 
The frequency dependence of g can be neglected for in- 
verse propagation times c/r between the two atoms much 
larger than the frequency differences lo' — u between the 
exchanged photons (which is the case for typical experi- 
mental parameters). Note that the integral over the fre- 
quency Lu of the photon emitted by the left atom expresses 
the correlation between both atoms induced through the 
photon exchange. 

Analogously, according to Eq. (3c), Pi{uj';ujd) rep- 
resents the emission of a negative frequency amplitude 
ujjy (incident on the detector) and a positive frequency 
amplitude lo (towards the other atom) by the left atom 
in Fig. [1]d, which, in turn, is subject to a negative fre- 
quency amplitude lo' originating from the right atom. 
The right atom, subject to a positive frequency ampli- 
tude Lu from the left atom, radiates a negative frequency 
amplitude w' (towards the left atom) and a positive fre- 
quency amplitude ujj (into the detector), represented by 
Pi{uj; LU]:)) in Eq. (3b). Energy conservation further en- 
forces Lo' = lul + llJo — uj. For the crossed contribution 
to the total scattered intensity at the frequency lod, this 
entails 



/c(t^D) = |.g|' 



■J dL0Pl{uJL 



+ ujD-Lo;iOD)Pi{i^\ujD) ■ (5) 



Once again, both atomic dipoles are conditioned on each 
other through the frequency of the exchanged photon. 
Finally, the total intensities measured by a broadband 
detector are obtained by integrating over the detected 
photon frequency, i.e., Il,c = J dLUolL.ci^^D)- 

By numerical integration of Eqs. Q and ([5]), we have 
studied the behavior of the elastic and inelastic compo- 
nents of the ladder and crossed terms for different values 
of the laser-atom detuning 6 — lol — loq (with o^o the 
atomic resonance frequency) and of the Rabi frequency 

(which is proportional to the laser amplitude \£l\)- 
In particular, some examples of the normalized inelastic 
ladder and crossed spectra for several parameters of the 
driving field are presented in Fig. |2] At small Rabi fre- 
quencies (see Fig. |2]^a)), the background and interference 
spectra are qualitatively very similar to the analogous 
quantities in the helicity preserving polarization chan- 
nel [lOJ |13j, and can be deduced from the two-photon 
scattering amplitudes. When ^ 7 (Fig. |2jb)), the 
spectral features of the ladder and crossed terms can be 
interpreted in terms of scattering and self-interference of 
photons emitted by one atom on the dressed state levels 
of the other atom |13j . Note that negative values of the 
ladder term in certain spectral regions arises from ab- 
sorption dominating emission for doubly scattered pho- 
tons; the total contribution, including single scattering, 
is always positive. 
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Figure 2: (Color online) Inelastic double scattering spectra as 
obtained from the two-atom master equation (solid line: back- 
ground, dashed line: interference contribution), compared to 
the predictions of our new approach, Eqs. (4]5l, based on 
single-atom Bloch equations with bichromatic driving (open 
circles: background, filled circles: interference contribution). 
Laser-atom detuning S = —67, Rabi frequency Q — (a) O.I7 
and (b) IO7. Perfect agreement is observed, which highlights 
the exact equivalence of both approaches. 



For comparison with these results for the single-atom 
response under bichromatic driving, we calculate the 
same quantities ~ the ladder and crossed double scatter- 
ing spectra - using the familiar two-atom master equa- 
tion approach [TTJ [T^l [T3]. As already mentioned, this 
approach provides a general framework for treating in- 
tense laser field-atoms interactions including interatomic 
correlations induced by the exchange of photons, but is 
restricted to small numbers of scatterers. Both, Ij^iujo) 
and Ic{lo]j), derived from the two-atom master equa- 
tion, are found to coincide perfectly with the results of 
Eqs. ( |4|5| , for all tested sets of the driving field parame- 
ters, see Fig. [2] Motivated by this agreement between the 
two apparently very different approaches, we succeeded 
to prove their equivalence analytically - at least for the 
elastic component of the scattered light [17]. We expect 
that the proof can be extended to the inelastic compo- 
nent, as suggested by the good numerical agreement ob- 
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served in Fig. [2] Details of the proof will be given else- 
where. 

Let us stress again that the reduction of the two-atom 
problem to single-atom Bloch equations is possible only 
in the case of large distance k^r ^ 1 (and hence \g\ <C 1) 
between the atoms. In this regime, both atoms exchange 
only single photons, what, in turn, allows to describe the 
field acting on each atom by a classical field, see Eq. ([T|). 
The same argument remains valid for a multiple scatter- 
ing process with more than two atoms. Provided that the 
atoms are far away from each other, there are again only 
single photons exchanged between the individual atoms. 
Moreover, the possibility that two photons emitted by 
the same atom meet again at another atom before leav- 
ing the sample can be neglected. In other words, all 
photons incident on a particular atom originate from dif- 
ferent previous atomic scatterers, and are hence uncor- 
related with each other. Since the quantum properties 
of the electromagnetic field manifest themselves only in 
correlations between different photons, the field inside 
the disordered sample can therefore be described by a 
classical field. This allows, in principle, to generalize the 
diagrammatic multiple scattering theory valid for clas- 



sical nonlinear scatterers (5] to the atomic case. The 
main difficulty to be addressed in future work will be to 
find an efficient method to solve the single atom optical 
Bloch equation for a fluctuating, polychromatic field, rep- 
resenting the radiation emitted by many randomly placed 
atoms. 

In summary, we have shown that the spectrum of in- 
tense laser light scattered by two distant atoms can be 
reproduced by solving the quantum mechanical time evo- 
lution equations for a single atom in the presence of a 
weak probe field representing the light emitted by the 
other atom. The interaction between the atoms is fully 
taken into account by integrating over the frequencies 
of the exchanged photons. Since, in contrast to other 
standard methods in theoretical quantum optics, it is 
not necessary to perform calculations in the composite 
Hilbert space of both atoms, our method will exhibit very 
favourable scaling with increasing system size, and can 
be applied to a larger number of atoms. Thereby, it can 
serve as a starting point for the development of a theory 
of multiple scattering of intense laser light by optically 
thick samples of cold atoms. 
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